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Abstract. We have investigated the zero and finite temperature behaviors of the 
anisotropic antiferromagnetic Heisenberg XXZ spin-1/2 chain in the presence of a 
transverse magnetic field (h). The attention is concentrated on an interval of magnetic 
field between the factorizing field (hf) and the critical one (he). The model presents 
a spin-flop phase ioi < h < hf with an energy scale which is defined by the long 
range antiferromagnetic order while it undergoes an entanglement phase transition 
at h = hf. The entanglement estimators clearly show that the entanglement is lost 
exactly at h = hf which justifies different quantum correlations on both sides of the 
factorizing field. As a consequence of zero entanglement (at h = hf) the ground state 
is known exactly as a product of single particle states which is the starting point for 
initiating a spin wave theory. The linear spin wave theory is implemented to obtain 
the specific heat and thermal entanglement of the model in the interested region. A 
double peak structure is found in the specific heat around h = hf which manifests the 
existence of two energy scales in the system as a result of two competing orders before 
the critical point. These results are confirmed by the low temperature Lanczos data 
which we have computed. 
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1. Introduction 

The zero temperature phase diagram (i.e. quantum phase diagram) of a model 
gives important information on the low temperature behaviors of the system [H |2]. 
The anisotropic antiferromagnetic Heisenberg (XXZ) spin-1/2 model shows different 
quantum phases with respect to a symmetry-breaking transverse field (non-commuting 
magnetic field) [31 HI [5]. The non-commuting field imposes quantum fluctuations into 
the ground state which can induce new phases. CS2C0CU is a quasi-one dimensional 
spin-1/2 XY-like antiferromagnet with weak inter-chain couplings [J'/ J = 0.014) which 
can be studied in terms of XXZ chain with anisotropy parameter A = 0.25 [6^, The 
scaling behavior and quantum phase diagram of the XXZ model in the presence of a 
transverse magnetic field (h^) have been investigated [3llll|5l|8]. 

Recently, the quantum spin models received much attentions from the quantum 
information point of views. These are prototype models to implement and examine 
the idea of quantum computations which requires the quantum correlations measured 
by entanglement. Hence, disentangled ground states have to be avoided for such 
implementations. On the other hand, the zero entanglement property of a ground 
state provides an exact form for it in terms of product of single particle states. The 
investigation in this direction resumed recently following the seminal work of J. Kurman 
and his collaboratores|9j. Several efforts have been devoted to this direction which is 
important for condensed matter researchers, i.e finding an exact (factorized) ground 
state even at particular values of the coupling constants [HI [121 113]. The factorized 
(exact) ground state is an accurate starting point to investigate the quantum nature of 
a phase close to the factorizing point in addition to some exact knowledge which gives at 
the factorized point. This property is implemented in this article to initiate a spin wave 
theory to describe the thermodynamic properties of the XXZ model in the presence of 
a transverse magnetic field. 

The U(l) symmetry of the XXZ model is lost upon adding the transverse magnetic 
field. Initially, a perpendicular antiferromagnetic order is stabilized by promoting a 
spin- flop phase (which has a partial moment projection along the field direction). At 
the factorizing field (hj) in the spin-flop phase the ground state is known exactly as a 
direct product of single spin states and the staggered magnetization along y-direction 
is close to its maximum value. In our model the factorizing field is hj = J^^/2{l^^~Kj, 
where quantum fluctuations are uncorrected and the ground state is the classical one (J 
is the scale of energy and A is the anisotropy parameter). For slightly larger magnetic 
field very close to the critical one {hf < h < he) the antiferromagnetic order becomes 
unstable and the staggered magnetization falls rapidly to vanish at the critical point. 
For h > he the spins become aligned in x direction and a fully polarized phase will be 
appeared. The region between the factorizing and the critical fields {hf < h < he) is 
the main issue of our study which is not well understood so far. 

The zero temperature properties of the intermediate region {hf < h < he) induces 
its signature into the thermodynamic functions of the model. We have found that 



Thermodynamics of XXZ chain around the factorizing field 



3 



this region is characterized by two energy scales and its fingerprint will appear as a 
double peak in the specific heat. Moreover, the existence of more than an energy scale 
in the model can be related to a spontaneous symmetry breaking (SSB) [TO]. Each 
broken phase is described by an order parameter which can be zero in the disordered 
phase. For the aforementioned model the symmetry breaking occurs where the staggered 
magnetization becomes zero. 

The thermodynamic properties of the spin 1/2 XXZ chain in a transverse magnetic 
field have been studied using the Low temperature Lanczos method [I5]. However, in 
this paper we will focus on the intermediate values of h and more precisely on the region 
where a double peak appears in the specific heat of the model. We implement the exact 
factorized ground state ai h = hf (where the entanglement vanishes) to build up a 
spin wave theory appropriate to describe the model in the intermediate values of the 
magnetic field. In the linear spin wave approximation we calculate the specific heat and 
thermal entanglement of the model. Moreover, the spin wave theory gives two different 
types of quasi particles which are representing the two energy scales. We have also 
studied both the zero and finite temperature properties of the model on a finite chain 
using the low temperature Lanczos method [14J. Our numerical results are in agreement 
with the spin wave theory counterparts. In the next section, we will briefiy review the 
zero temperature properties of the model from the quantum information point of view. 
We will provide the spin wave theory in Sec J3] where the quantum property, specific heat 
and thermal entanglement are obtained. Finally, the numerical Lanczos results will be 
presented in addition to discussions on the mentioned topics. 

2. Zero temperature phase diagram 

The anisotropic spin-1/2 Heisenberg model in the presence of a transverse magnetic field 
is described by the following Hamiltonian, 

H = jY.{s:s:^, + sfsf_,, + As^su, + hs:), (i) 

i 

where 5*°" s are spin-1/2 operators, A is the anisotropy parameter, h is proportional to 
the transverse magnetic field and the exchange antiferromagnetic coupling J which 
defines the scale of energy is set to one. This model has several quantum phases 
with respect to the transverse magnetic field and anisotropy parameter (A). The a- 
component spin structure factor at momentum p is defined by 

N 

G-»{p) = Y.iS^S'^^Je^P^, (2) 

x=l 

and its increasing behavior in term of system size shows the magnetic order at the 
specified momentum. We have plotted in Fig.([T]) the magnetization along x direction 
(Mj = (1/A^) ^('S'f )) and the y-component of spin structure factor at momentum p = n 
versus h. The value of A = 0.25 has been fixed to fit the case of real material CS2C0CU. 
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The numerical data has been obtained by zero temperature Lanczos method on a finite 
chain of length = 20 with periodic boundary condition. 
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Figure 1. (color online) Results of numerical zero temperature Lanczos method for 
the magnetization along a;-direction (right vertical axis) and the y-component spin 
structure factor at momentum vr (left vertical axis) versus the transverse magnetic 
field h for a chain with = 20 and A=0.25. The magnetization in z-direction is zero. 

The magnetization curve can be distinguished in three parts: (1) Magnetization 
for Q < h < hf, (2) Magnetization for hf < h < he, and (3) the paramagnetic phase 
for h > he- At zero magnetic field, there is no order in the system and all order 
parameters are zero i.e M^^y^z = 0. The nonzero value of h starts to align the spins in x 
direction and induces a small ferromagnetic order in x direction. The magnetization, M^, 
increases monotonically by increasing the magnetic field. For h < hf the quantum effects 
are considerable and the magnetization changes parabolically versus h. Increasing h, 
suppresses the quantum correlations and they goes to zero around hf. For hf < h < h^ 
the magnetization is increased linearly versus h. Magnetization increases up to the 
critical point (he) and saturates at infinite field. Although the staggered magnetization 
along y-direction becomes zero at the finite critical field, the magnetization in x-direction 
will be fully saturated only for the isotropic case (A = 1) [13]. In other words, the full 
saturation in x-direction for A 7^ 1 will happen for h ^ 00. For h > h^ all of spins 
align almost completely (for A 7^ 1) in the x direction and we have a fully polarized 
paramagnetic phase. 

We would like to draw your attentions to the region hf < h < h^, where the model 
behaves surprisingly. Let us first study the spin chain through the entanglement {T2) of 
two spins [ini [T7] which have no classical counterpart. The entanglement is defined by 
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the following relation, 

Where, Cij is the concurrence [H] which is implemented instead of the pairwise zero 
temperature entanglement between two spins at sites i and j. For XXZ model, is 
zero thus the concurrence takes the form [TTj, Cij = 2 max{0, C^j\ C'^-J^}, where 

cS^ = {s:sj) + \{sfs])-{S!s^)\-\, 

cf^ = Ksfs]) + {s^s;)\ - ^J{\ + {SfSpy-{M^r. 

(4) 

Using quantum Monte Carlo simulation, T. Roscilde and collaborators [10] have 
shown that unlike the standard magnetic order parameters (FigUI) the pairwise 
entanglement, plays an essential role at the factorizing field. At the factorizing field 
(also called classical field, hd = hf = \/2{l + A)) [9] the ground state takes a 
product form [9l [13] and its entanglement is zero. In Fig.(|2]), we have plotted Cj''^^^ 

(2') 

and Cl/_^i versus the transverse field h for XXZ spin-1/2 chain with A = 0.25 by 
using exact diagonalization Lanczos method. At the classical field {h = hf c:^ 1.58), 
Ci^i+i = Cf^^]^ = C^'f_^i = 0. At this point the ground state of the model is factorized 
(disentangled) where the Neel order is roughly maximized in y direction. 
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Figure 2. (color online) Entanglement estimators of the XXZ spin-1/2 chain versus 
transverse field h with parameter A = 0.25. C*^^-' is depicted by (black) circles and 
C(2) by (red) squares. At the classical field {h = hf ~ 1.58), c[%^ = 0*^^^ = 0. 

Let us describe the behaviors of C^j^ and C^j^ from a spontaneous symmetry 
breaking (SSB) point of view. It is found in Ref. [lO] that the competition between the 
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two functions C^p and CfJ"* demonstrates the appearance of a SSB. In other words, 
for our model, a SSB can be occurred when CfJ^ < CfP. At this condition the 
magnetic field is greater than the factorizing value {h > hf). Moreover, SSB is usually 
recognized to happen at the position where the order parameter becomes zero. In our 
model two standard order parameters, and SMy (staggered magnetization along y 
direction) can represent the quantum phases of the model. We have plotted in Fig.(Il]) 
the magnetization along x direction (M^.) and the y-component of spin structure factor 
(G^^) at momentum vr for A = 0.25 XXZ spin-1/2 chain. is nonzero in the whole 
range of magnetic fields however, G^'^{tt) which is used to show the antiferromagnetic 
order of the system has a different behavior. As it is observed from our Lanczos data, 
G™(7r) increases hj h up to hf. This behavior is also found in the SM'^ curve which 
has been obtained by density matrix renormalization group (DMRG) [8]. Thus, no SSB 
occurs for h < hf. By increasing the magnetic field for h > hf, G^^(7r) (or equivalently 
SMy) decreases rapidly and falls to zero at the critical point h^. In this respect, a 
symmetry breaking can occur only for h > hf. In this region, as mentioned before the 
slop of with respect to h is different from the corresponding one for h < hf. 



3. Spin waves theory 

As discussed in the previous section the entanglement is zero at the factorizing field 
(hf). It allow s us to write the many body ground state as the direct product of the 
single spin states on odd and even sublattices |9] 

\GS)= (g) \S,)\S,) (5) 

i£odd,j^even 

The spin state on each sublattice is expressed in terms of polar angles {6, (p) which 
defines the rotation of the up-spin eigen-state of to the specific direction defined 
by the factorized state. Let us label the odd (even) sublattice by A{B). Thus, (6*, 0) 
represents A-sublattice while (/3, a) is the corresponding one for B-sublattice. It has 
been shown [13j that one can consider = = 0; without loss of generality. Moreover, 
the magnitude of remaining polar angles are equal \6\ = \/3\ in the case of homogeneous 
spin model (like here) and is given by 

cos(^) = -/^ (6) 

For the antiferromagnetic model the factorized ground state is defined by /3 = —6. 

Before starting the spin wave approach we implement a unitary transformation 
on the Hamiltonian. All spins on the A-sublattice are rotated with angle 6 
counterclockwise around y-direction and clockwise for spins on B-sublattice. The rotated 
Hamiltonian (H) is the result of rotations on all lattice points, H = D^HD and 
D = <S)ieAjeB The single spin rotation operator is T>i{9) = exp{—i9Sf /h). 

In the rotated spin model, the Hamiltonian will have a factorized ground state which is 
a ferromagnet at the factorizing field. The rotated spin Hamiltonian is written in terms 
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of boson operators a, b with the following Holstein-Primakoff transformation, 

Sj^j^ = (215 ojctj) ^ fli, Sj^j^ = S 

where Sa{b) = 1^a(^b)^MB)'^a(b) are the rotated spin operators and T>a(b) is the 
unitary single spin rotation operator. The bosonic Hamiltonian in the linear spin wave 
approximation, i.e. Sj^^ ~ y/2Sai; S^^ ~ V2Sbi, is given by 



Ah ^ h 

- f + + h.c] + 1(A _ l)v/r^(a, + h + h.c)). (7) 



To diagonalize the bosonic model we first implement the Fourier transformation and 
then apply a rotation to the boson operators (a^, 6^) 



6, = ^5^e-('='+t)6, 



il^k = cos rjkak - sin r]kbk, 

Xk = smT]kak + cos?7fc6fc. (8) 

The diagonalized Hamiltonian in terms of two sets of quasi-particle operators is given 
by 

n = Eo + Y^iujtxlxk + uJk^i^k), (9) 

k 

where the excitation spectrums have the following forms 
a;± = -^(l + A)-A±Acos|, 
N 



E, = -{A- hhf) + uUt^-f + v/2iV(l-A2)(- - l)t' 



'f(l-A2)(l- A) 
t+ = ^ ^. (10) 

^0 

In the above calculations, a translation Xo ^ Xo + has been performed in the 
diagonalization procedure of the Hamiltonian ((Tj). 

The Hamiltonian is now represented in terms of two quasi-particles (bosons), each 
defines an energy scale. The two energy scales which are the excitation energies of each 
bosons lead to two different dynamics for the system. The consequence of two types of 
dynamics will be seen in the structure of specific heat which will be studied. However, 
the existence of two different quasi-particle energies (cu^ ^k) ^^^^ ^^S^ 

dynamics in the model. 
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3.0.1. Specific heat To get the finite temperature properties of tlie model, we assume 
(n^,n ), for tlie spin wave distribution functions, where Ph(n~^,n ) 



is the probabihty of parallel {n^ = xlxk) and perpendicular (n^ = iplipk) normal modes 
appearing in the fc-momentum state which satisfies Yln+ n- Pk{n^ ,fi~) = 1 for all fc's. 
The substitutions of = (xlxk) and fi^ = (iplipk) (where (■ ■ ■) represents the thermal 
average) in the spin-wave Hamiltonian ([9]) gives the free energy, 

F = Eo + J2i^M + ^kf^k)+Tj2 E Pk{n-,n+)\nPk{n-,n+). 

k k n+,n^ 

The number of bosons are controlled by the following constraint which is the 
magnetization in x-direction, 

^^ = ^-^E(-^+-^-)-^- (11) 

k 

The free energy is minimized with respect to Pk{n~^,n~)s under the constraint of (fTTi) 
which is applied by a Lagrange multiplier (/i) via the boson's occupation number 

= 1,2, • (12) 

The constraint flTT]) is applied by the values of Mrc{h,T) which have been obtained by 
the numerical Lanczos method. We have plotted in Fig. ([3]) the specific heat of A = 0.25 
XXZ spin-1/2 chain versus temperature and different values of magnetic field. A double 
peak is observed in the specific heat which is the signature of the existence of two 
comparable energy scales. More precisely, the specific heat for h = 1.4 and 1.5 have 
a narrow peak at low temperature (T ~ 0.2) and a broaden one at higher T. We will 
discuss more on this point in the next section. 



3.0.2. Thermal entanglement The established spin wave theory close to the factorizing 
point is used to obtain the thermal behaviors of the correlation functions in Eqs.(|l]) 
and consequently to calculate the concurrence. This method can describe the thermal 
entanglement of two spins in linear spin wave approximation. In this approach, one can 
find the following expression for Clf and C^]^ 

k 

k 

As it is seen from these functions, C^j^ is always less than zero thus the concurrence is 
2 max{0,C-j }. We have plotted in Fig.(jl]) the thermal entanglement of the A = 0.25 
XXZ spin-1/2 chain versus transverse field at temperature T = 0.05. By increasing 
the magnetic field, quantum fluctuations are decreased and the thermal entanglement 
is declined. At the factorizing point, quantum fluctuations become approximately 




Figure 3. (color online) The spin wave results for specific heat of the XXZ model 
versus T for different values of the transverse magnetic field h and anisotropy parameter 
A — 0.25. (Black) circle is for h = 1.2, (red) square: h = 1.4, (green) gradient: h — 1.5 
and (blue) diamond: h = 1.6. A narrow peak at low temperature followed by a broaden 
one for higher temperature are observed for h = 1.4, 1.5. 



uncorrelated and the thermal entanglement is very close to zero (~ 10"'^). It means 
that the thermal entanglement is mainly originated from the ground state and the 
excited states have very tiny contribution to the entanglement. 




h 

Figure 4. Thermal entanglement of the XXZ model versus h and anisotropy parameter 
A = 0.25 at T = 0.05 which has been obtained by spin wave approximation. Because 
of the tiny contribution of excited states to the entanglement, its value at /i/ is ~ 10~^. 
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We have studied the effects of transverse magnetic field on the zero and finite 
temperature properties of A = 0.25 XXZ spin-1/2 chain. We have focused our attentions 
on the intermediate region oi hj < h < he where the model behaves more interestingly. 

M. Kenzelmann and his collaborators [6] have investigated experimentally 
the effects of a transverse magnetic field on the quasi-one dimensional spin-1/2 
antiferromagnet CS2C0CU, using single- crystal neutron diffraction. Due to the weak 
inter-chain couplings in CS2C0C/4 {J'/ J = 0.014) [7] where J is the coupling within a 
chain, it is proposed in Ref. [6] that the bulk material has a spin liquid phase at the 
interval hf < h < h^. In a spin liquid phase all order parameters should be zero and 
there should be no long range order in the system. The existence of very weak coupling 
between magnetic chains in this material makes it feasible to be described by a one 
dimensional spin-1/2 XXZ model with the anisotropy parameter A = 0.25. However, 
for the ID XXZ model in the region hj < h < he both magnetization and staggered 
magnetization are nonzero. Thus the ground state of the model could not be a spin 
liquid phase. Although the behaviors of the magnetic orders for this region are clearly 
known, the lacuna of a perfect study of the properties of the system at the intermediate 
region is still felt. In this respect, we have devoted our attentions to survey the magnetic 
and thermodynamic behavior of the system at the intermediate region of the transverse 
field. 

We have implemented the low temperature Lanczos method (LTLM) [TJl [15] to 
compute the thermodynamic behaviors of the model for a chain of finite length. LTLM 
has been used since it is accurate for low temperatures and specially the thermodynamic 
averages reach the ground state expectation values as temperature approaches zero. The 
specific heat versus temperature for different values of the magnetic field on a chain of 
= 20 and with A = 0.25 has been plotted in Fig. ([5]). The sampling is taken over 
R = 100 Lanczos numerical data which have been obtained by different initial random 
states. The finite size effect can be ignored since the number of sampling (R) is rather 
high. The result of the spin wave theory (FigjS]) and LTLM (FigjS]) are in mutual 
agreement; moreover, both figures show the presence of double peak in the specific 
heat for h = 1.4, 1.5 which is an evidence for the existence of two scales of energy or 
equivalently two dynamics in the system. 

The specific heat versus T for all values of the magnetic fields has Schottky like peak 
at low temperatures which is a remarkable feature of the antiferromagnetic behavior. 
For small values of the magnetic field {h < 1.4) the Schottky anomaly is the only one 
which justifies the existence of a single dynamics in the model. Further increasing of 
the magnetic field (being close to the factorizing point) a broaden peak is emerged in 
the specific heat data versus T. This bump is the result of the paramagnetic order. 
Two different orders are in competition and its signature is specially observable for 
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Figure 5. The numerical Lanczos results for the specific heat of XXZ model versus 
T for different values of transverse magnetic field and A = 0.25. At the intermediate 
region of h, a double peak structure is observed which indicates the existence of two 
energy scales in the system. 

magnetic fields in the region hf < h < he where the magnitude of the two types of 
ordering becomes comparable. The antiferromagnetic order in y-direction is the result of 
exchange coupling in the broken U(l) symmetry phase. The effect of transverse magnetic 
field as a paramagnetic order in x-direction shows its presence when it is enough large 
to define a new scale of energy. Around the factorizing field the two orders manifest 
their influence on the model as a narrow and broaden peak of the specific heat which 
happens for hf < h < h^. This is the region where the spontaneous symmetry breaking 
is started to happen. At the critical point (he) the antiferromagnetic order vanishes 
and paramagnetic order is the only representative of the model. A broaden peak in the 
specific heat versus T is significant for h > he which justifies the paramagnetic order. 

Employing the specific heat data, one can also scan the behavior of the energy gap 
(Eg). The existence of energy gap in the model is clearly observed by the exponential 
decay of the specific heat at very low temperatures (T < 0.02). At enough low 



temperatures (T < Eg) the specific heat and energy gap are related by [7] 

c^T^/^ oc e""^. (13) 

Thus, the slop of (c^T^/^) curves versus ^ in log scale gives good information about 
the energy gap. In Figs.( l6|7|) we have plotted {c^T^^'^) versus i for different values of 
transverse field. Fig.(l6]) shows the gap behavior close to the factorizing point, h < hj. 
While its behavior for h > hf is presented in Fig. ([7]). A turning point of different plots 
in Fig. ([7]) for different transverse magnetic field is the result of different behavior for 
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Figure 7. c^T'^/^ of XXZ model versus l/T for different values oih > hf and A = 0.25. 
At low temperatures, the slop of the curve is decreased from h = 1.45 to /i = 1.6 and 
then starts to increase by increasing of ti. 
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h < he — 1.Q5 (decreasing gap) and h > he (linear increasing paramagnetic gap). Our 
results fit very well with Figs. (5.14 and 5.15) of Ref. [7]. 

Motivated by neutron-scattering results [6], T. Radu [7] investigated the effects 
of non-commuting field on the ground state of CS2C0CI4. To be more precise in 
comparison, let us refer also to Fig. 13 of chapter. 5 of Ref. [7] where the experimental 
data of the specific heat have been shown. The experimental data (specially in the 
inset of Fig. 13 of chapter. 5 of Ref. 0) display a double peak in the specific heat versus 
temperature which is for the intermediate range of magnetic field (less than the critical 
one). 

It is also worth to point out the behaviors of the system from the internal energy 
points of view. The internal energy and specific heat are related by relation C = The 
double peak structure of the specific heat presages that at the intermediate values of the 
transverse field, there are two temperatures where the internal energy of the system gets 
its maximum variation. In other words, at these temperatures the maximum amount of 
states contribute to the response functions of the systems. Thus, one can also conclude 
that the appearance of the two energy scales in the system is appropriate with the 
number of contributed states. Knowledge on these properties could be important in the 
study of the magneto-caloric effects of the XXZ model in the transverse field which is a 
work in progress. 
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